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Abstract. In [16) Tiwari proved that termination of linear programs 
(loops with linear loop conditions and updates) over the reals is decid- 
able through Jordan forms and eigenvectors computation. In [3] Braver- 
man proved that it is also decidable over the integers. In this paper, 
we consider the termination of loops with polynomial loop conditions 
and linear updates over the reals and integers. First, we prove that the 
termination of such loops over the integers is undecidable. Second, with 
an assumption, we provide an algorithm to decide the termination of a 
class of such programs over the reals. Our method is similar to that of 
Tiwari in spirit but uses different techniques. Finally, we conjecture that 
the termination of linear programs with polynomial loop conditions over 
the reals is undecidable in general by reducing the problem to another 
decision problem related to number theory and ergodic theory, which we 
guess undecidable. 



1 Introduction 

Termination analysis is an important aspect of program verification. Guaranteed 
termination of program loops is necessary for many applications, especially those 
for which unexpected behavior can be catastrophic. For a generic loop 

while (conditions) {updates}, 

it is well known that the termination problem is undecidable in general, even for 
a simple class of polynomial programs [3]. In [5] Blondel et al. proved that, even 
when all the conditions and updates are given as piecewise linear functions, the 
termination of the loop remains undecidable. 

In [16] Tiwari proved that termination of the following programs is decidable 
over R (the real numbers) 

P : while (BX > b) {X := AX}, 

where A and B are respectively n x n and m x n matrices, BX > b represents 
a conjunction of linear inequalities over the state variables X and X :— AX 
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represents a (deterministic) linear update of each variable. Subsequently in [3] 
Braverman proved that the termination of P is decidable over Z (the integers) . 
In this paper, we consider the problem of termination of the following loop: 

Pi : while (P(X) > 0) {X := AX}, 

where X = [x\ ... xn] t is the vector of state variables of the program, P{X) = 
[Pi(X) P2PO • • • Pm(X)] T are polynomial constraints, each Pi{X) (1 < i < m) 
is a polynomial in Q[X] and A is an N x TV matrix over Q (the rational numbers). 
That is to say, we replace the linear constraints in P with polynomial constraints 
and keep linear updates unchanged. 

There are some well known techniques for deciding termination of some spe- 
cial kinds of programs. Ranking functions are most often used for this purpose. 
A ranking function for a loop maps the values of the loop variables to a well- 
founded domain; further, the values of the map decrease on each iteration. A 
linear ranking function is a ranking function that is a linear combination of 
the loop variables and constants. Recently, the synthesis of ranking functions 
draws increasing attention, and some heuristics concerning how to automati- 
cally generate linear ranking functions for linear programs have been proposed 
in |8l9ll2j . In [T2] Podelski et al. provided an efficient and complete synthe- 
sis method based on linear programming to construct linear ranking functions. 
In [5] Chen et al. proposed a method to generate non-linear ranking functions 
based on semi-algebraic system solving. However, existence of ranking function 
is only a sufficient condition on the termination of a program. It is not difficult 
to construct programs that terminate, but do not have ranking functions. 

To solve the problem of termination of Pi, we do not use the technique of 
ranking functions. Our method is similar to that of Tiwari in spirit. Our main 
contributions in this paper are as follows. First, we prove that the termination of 
Pi over Z is undecidable. Then it is easy to prove that, if ">" is replaced with ">" 
in Pi, termination of the resulted Pi over Z is undecidable either. Second, with 
an assumption, we provide an algorithm to decide the termination of Pi over KL 
Finally, we conjecture that the termination of Pi over R is undecidable in general 
by constructing a loop and reducing the problem to another decision problem 
related to number theory and ergodic theory, which we guess undecidable. 

The rest of the paper is organized as follows. Section 2 proves the undecid- 
ability of Pi and its variation over Z. Section 3 introduces our main algorithm. 
The main steps of the algorithm are outlined first and some details of the steps 
are introduced separately in several subsections. With an assumption, we prove 
the correctness of our algorithm at the end of Section 3. After presenting our 
conjecture that the termination of Pi is generally undecidable in Section 4, we 
conclude the paper in Section 5. 

2 Undecidability of P t over Z 

Definition 1. A loop with N variables is called terminating over a ring R if for 
all the inputs X £ R N , it is terminating; otherwise it is called nonterminating. 



The undecidability of Pi is obtained by a reduction to Hilbert's 10 th problem. 
Consider the following loop: 

P 2 : while {x m+1 - f{xi,... 1 x m f >0) {X := AX} 

where X ~ \x\ ...x m+ i] T , A = diag(l, 1, 1/2) is a diagonal matrix and 
f(x%, . . . , x m ) is a polynomial with integer coefficients. 

Lemma 1. For any input (x%, . . . , x m +\) € Z m+1 , P 2 terminates if and only if 
f(xi, . . . , x m ) does not have integer roots. 

Proof. (=>) If f(xi, . . . , x m ) has an integer root, say (yi, . . . , y m ), obviously Pi does 
not terminate with the input Y = (yx, ■ ■ ■ , y m , 1)- 

(•<=) If f(xi, . . . ,x m ) has no integer roots, for any given X G Z m+1 , —f(xi, . . . , x m ) 2 
is a fixed negative number. Because (xi, x m ) will never be changed and (l/2) n — > 
as n —* +oo, the loop will terminates after sufficiently large n iterations. 

Theorem 1. Termination ofP± over "L is undecidable. 

Proof. Because the existence of an integer root of an arbitrary Diophantine 
equation is undecidable, the termination of Pi over Z is undecidable according 
to Lemma [TJ 

If ">" is substituted with ">" in Pi, the loop becomes 
Pi : while (P(X) > 0) {X := AX}. 

Analogously, we denote by Pj the loop obtained by substituting ">" for ">" in 
P2. It is easy to see that Lemma [T] still holds for P' 2 . Then we get the following 
theorem. 

Theorem 2. Termination ofP[ over Z is undecidable. 

3 Relatively Complete Algorithm for Termination of Pi 
over M 

To decide whether Pi terminates on a given input X £ M. N , it is natural to 
consider a general expression of A n X, for instance, a unified formula expressing 
A n X for any n. If one has such unified formula of A n X, one can express the values 
of P(X) after n iterations. Then, for each element of P(X) (each constraint), i.e., 
Pi(X), one may try to determine whether Pi(X) > as n — > +00 by guessing 
the dominant term of P%(X) and deciding the sign of the term. 

That is the main idea of our algorithm which will be described formally in 
Subsection 13.21 At several main steps of our algorithm, a few techniques and 
results in number theory and ergodic theory are needed. For the sake of clarity, 
the details are introduced subsequently in the next subsections. 

The first subsection is devoted to expressing A n X in a unified formula. 



3.1 General Expression for A n X by Generating Function 

In this subsection the general expression of A n X will be deduced with generating 
function, not Jordan form. 

Lemma 2. |13j Let ot\, . . . , ad be a sequence of complex numbers, d > 1 and 
ad 0. The following conditions on a function f : N — > C are equivalent to each 
other: 

i. ^2 f{ n ) xn — qQ , where Q(x) — l+aix+. . .+adX d and P{x) is a polynomial 

n>0 

in x of degree less than d. 

ii. For all n>0, f(n + d) + gji/(b + d - 1) + a 2 f(n + d-2) + ... + a d f{n) = 0. 

k 

iii For all n > 0, f(n) — ^ Pj(n)7™, and Q(x) = 1 + a\x + a 2 x 2 + . . . + adx d = 

i=i 

k 

EI (1 — lix) di , where the 7$ 's are distinct, and Pi{n) is a polynomial in n of 

i=l 

degree less than di. 

Corollary 1. A is a d x d square matrix with its entries in Q. Suppose that 
the characteristic polynomial of A is D(x) = x d + aix d ~ 1 + . . . + a,;_ s i s , where 
<*d-s 7^ and s > 0. Define f(n) — A n+s X and let fj(n) be the j-th component 
of f(n). Then for each j, fj(n) can be expressed as 

k 

/j(n) = EwW{T 1 (1) 

i=l 

where £j 's are all the distinct nonzero complex eigenvalues of A and pji (n) is a 
polynomial in n of degree less than the multiplicity of £i . 

Proof. First, A d + a\A d ~ x + ■ • ■ + a d _ s A 3 = because D(x) is the characteristic 
polynomial of A. So, for any n > 0, 

f(n + d-s) + axf{n + d - (s + 1)) + • • • + ot d ^ s f(n) 
= A n+d X + a 1 A n+d_1 X + ■ ■ ■ + a d - a A n+s X 
= (A d + a^- 1 +■■■ + a d - s A 3 )A n X = 0. 

Thus, for each j, fj (n + d — s) + ai/ 3 (n + d — (s + 1)) H + a d - s fj (n) = 0. By 

fe fe 
LemmaH / f (n) = £ Pii(™)£" and Q(x) = 1 + am + h ct d - 3 x d - s = FJ (1 - £ix) d * 

i=l i=l 

where Pjiin) is a polynomial in n of degree less than di. It's easy to see that x = 
is not a solution of Q(x) and J3* =1 di = d~ s. Because 

D(x) = = x d ~ s + a^-*' 1 + ... + a d _ s = x d ]J(l - =x s l[(x- &) d « , 

i=l '' i=l 

£i's are all the distinct nonzero complex eigenvalues of A and di is the multiplicity 
of £i. That completes the proof. 



Remark 1. According to Corollary [TJ we may compute the general expression 
of A n X as follows. First, compute all the complex eigenvalues of A and their 
multiplicities. Second, suppose each fj(n) of f(n) is in the form of eq.([T]) where 
the coefficients of pji are to be computed. Third, compute /(l), . . . , /(d), and 
obtain a set of linear equations by comparing the coefficients of the resulted 
fi(i) (1 < i < d) to those of eq.([TJ). Finally, by solving those linear equations, we 
can obtain fj(n) and /(n). 

Example 1. Let's consider the following loop: 

while (x 5 + x\ + x\X2 — x\ — 2x3X4 — x\ > 0) {X := AX; }, 

where 
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We shall show how to compute the general expression of A n X by taking use of 
Corollary [TJ The characteristic polynomial of A is 

D(X) = A 5 + —A 4 + —A 3 + ^A 2 + — A + - = (A + -)(X 2 --X + 1)(A 2 + A + 1). 
VJ 10 10 5 10 2 v 2 A 5 A ; 

The eigenvalues of A are 
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Set f(n) = A n X = [A(n) / 2 (n) / 3 (n) / 4 (n) / 5 (n)] T . For any n > 0, we have 
/(" + 5) + ^/(n + 4) + Lf{n + 3) + + 2) + JL/(n + 1) + \f(n) = 0. 

Because the multiplicities of £1,62, £3, £4 and £ 5 are all 1, by Corollary [TJ we may 
assume for 1 < j < 5 

55555 
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tively, and by solving some linear equations (see Remark [JJ we can obtain 
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2),/(3),/(4) and /(5) equal to AX, A 2 X, A 3 X, A 4 X and A 5 X respec- 
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3.2 Main Algorithm 



According to subsection l3.1i the general expression of A n+m X is a polynomial in 
xn,u and £i> • • • > £g> the nonzero complex roots of D{x). If we substitute 
A n+m X for X in P(X) and denote the resulted Pj(X)(l < j < m) £ P(X) by 
Pj (X, n) , then Pj (X, n) can be written as 

Pj{X,n) = Pj o(X,n)+ Pjl (X,n)rii + ■ ■ ■ + p 3M (X,n)rfi A , (2) 

where rjk(l < k < M) is the product of some £j's. 

To determine whether Pi terminates, we have to determine for each j whether 
Pj(X,n) > holds for all n. To this end, it is sufficient to know whether the 
dominant term (leading term) of eq.([3]) is positive or not as n — > +oo. In the 
following, we shall give a more detailed description of eq.@ so that we can 
obtain the expression of the leading term of eq. © . 

Let r\k — rfce"* 2 *" 1 , where i = T and r^ is the modulus of rjk- Without loss 
of generality, we assume r\ < T2 < ■ ■ ■ < tm- For convenience, set T]o = To = 1. 
Rewrite Pj (X, n) as 

Pj(X,n) =pjo{X,n)r + pji(X,n)e r 1 H h p jM {X,n)e M r M . 

Suppose Tj is the common period of all the e aq27T ' (1 < q < M) where a q is a 
rational number. 

m s — 1 

Definition 2. For each j (1 < j < £ T ; ), i/ j = E T « + * a™d 1 < » < T s , then 

i=i i=i 

define 

Gj(X,n) = P s (X,T s n + i-l). 

m 

Notation 1 For each i (1 <i < Ti), expand Gj(X,n), collect the result with 

i=i 

respect to (w.r.t.) n l r^, and let Cjki(X,n) denote the coefficient of the term n l r^. 
Then Gj (X, n) can be written as 

C jw (X,n)r? + Cju(X,n)nr? + ■■■ + C Jldl (X,n)n dl r? + ■■■ (3) 
+CjMo(X,n)rM + CjMi(X,n)nrM H + C jM d M {X,n)n dM >m, 

where di(l < I < M) is the greatest degree of n in Gj(X,n) w.r.t. ri. 

It can be deduced that if < rj, the order of n ll r™ is less than the order of 
n' 2 r™ for any l\ and 1% when n goes to infinity. Similarly, if li < h, the order of 
n ll r™ is less than the order of n l2 r™. So, it is natural to introduce an ordering 
on the terms n l r™ as follows. 

Definition 3. We define n ll r" ~< n l2 r r f if ri < rj or r*j = rj and l\ < l^. A term 
Cjkin 1 ^ in eq. ([3]) is said to be the leading term and Cjki the leading coefficient 
ifn l r"j occurring in Cju is the largest one under that ordering -<. 



Suppose Gj(X,n) = P s (X,T s n + t) for some s and t. For those e aq2wi, s where 
Qq's are rational numbers, e ( T s"+ t ) a <i 2 ' rl — g^ 2 ™ because T 8 is the common pe- 
riod. Because there may be some e ( T = n + t )° ! i 2nl, s with irrational a q 's, each Cjki(X, n) 
can be divided into three parts, 

Cjki(X, n) = Cjkio(X) + Cjkii(X) + Cjki2(X, n), 

where C jk io(X) does not contain any e ( T ^+^ q ^\ C jk n(X) contains e (Tsn+t)aq2ni 
with rational a q and Cj k n{X,n) contains those e ( T = n + t ) a q 2 ^ 1 with irrational aJ3- 
Further, C jk i2(X,n) can be written as 

C ]k n(X,n) = C ik ii(X, sin((nT a + i)a fc i27r), cos((nT s + t)a fc i27r), . . . , 
sm((nT s + t)a fcSfc 27r), cos((nT s + t)a kSk 2-k)), 

where {a k i2-K, . . . , Qfe Sfc 27r} is a maximum rationally independent groupH. 

Example 2. We continue to use the loop in Example 1 to illustrate the above 
concepts and notations. 

Because |£i| = |f 2 | = 1, let £i = e ai27ri and £2 = e~ ai2ni , where a^n is the 
argument of £1. It's not difficult to check that on is an irrational number|f| 
For the sake of clarity, in the following we firstly reduce the expressions of 
/i( n )i fzin), h(ri), fi{n) and /s(n), and then substitute them in the loop guard. 
Let a 2 = |, then £3 = e a ^\ £4 = e"" 2 ^ 1 , and fi(n), / 2 (n), f 3 (n), h{n) and / s (n) 
can be rewritten as 

h(n 
h{n 

Mn 

U{n 



xi cos(nQi27r) + Xl X2 sin(7iai27r), 
X2 cos(nai27r) + - ^ sin(nai27r), 

X3 cos(nQ2i7r) H sm(na227r), 

, . \/3a;3 + v / 3s4 . , „ s 
X4 cos(na22n) sm(na?227r), 



Substituting /i(n), fa{n), fs(n), /4(n) and /s(n) for xi, x 2 , 2:3, 2:4 and 2:5 respec- 
tively in the loop guard, we get that the resulted loop guard is 

Lir^ + (L 21 +L 22 + L 23 )r^, 

where n = \,r 2 = 1, and 
Li = (-l)"x 5 , 

x\ + 2x3X4 + Ax\ 5a;i + x 2 — 2x1X2 



L21 = 



2 4 

„2 



2X4 — 2X3X4 — Xo 1— 1— 2 X , 

L22 = cos(n2a227r) — (V0X3X4 + V 02:4) sm(n2a227r), 

x? + xl - 6x1x2 , . 7xf - x\ - 2xix 2 . , . 
-L23 = cos(niai27r) H sm(n2ai27r). 



1 Later it will be proved that Cj k io{X), Cj k n(X) and Cj k i2{X,n) are reals for any n. 

2 "Rationally independent" will be described later. 

3 A general algorithm for checking whether an argument is a rational multiple of ir 
will be stated in detail in subsection 3.4. 



Since a 2 = §, the period of £| = e 2a227T = cos(2q 2 2tt) + sin(2a 2 27r)i is 3. Then we 
compute 

Gi(X,n) = ((-l) 3 "x 5 )r 3n + Ciaori", 
G 2 (X,n) = ((-l) 3n+1 x 5 )r 3 " +1 + C 220 r 3ll+1 , 
G 3 (X,n) = ((-l) 3 " +2 x 5 )r?" +2 + C 320 r 2 3n+2 . 
Take Gi (X, n) as an example. 
G120 = G1200 + G1201 + G1202, 

£3 + 2x3X4 + 4^4 5xf + X? 2 — 2X1X2 

C1200 = — ■ 1 - , 

G1201 = 



2? + xl - 6x1X2 . . 7x 2 -xl~2x 1 x 2 . , . 
Ci 2 o 2 = cos(Sn2ai2ir) H sin( 6n2ai2-K). 

Notation 2 We denote by Cj k i(X,n) y (and call Cj k i(X,n) "positive") if 

mm{C jk w(X) +Cjkii(X) +C'jki2(X, j/n, j/i 2 , . . . , J/ Sfc i, y Sfc2 )} > 

subject to {j/ 2 ! + y\ 2 = 1, • • • , y 2 fc i + y 2 fc2 = 1}. 

Remark 2. It is not difficult to see C jk i(X,n) y iff 

V(yn,yi2, • • ■ ,y, k i,y 3k 2) ((2/11 + 2/12 = 1) A • • • A (y 2 fcl + 2/ 2 fe2 = 1) -> 
G jfc!0 (X) + G 3fcn (X) + C jk i2{X,yu,yi2, ■ ■ ■ , y„ k i, 2/ Sfc2 ) > 0) 

because {(2/11,2/12, • • • ,y Sfc i, J/ Sfc2 ) : y?i + yf 2 = h i = 1> ■■■> is a bounded closed 
set. If y and > are replaced with y and > respectively in the notation, we get 
the notation of C jk i(X, n) y ("nonnegative"). 

Roughly speaking, for any Gj(X,n), if its leading coefficient C jk i{X,n) y 0, 
there exists an integer Nj such that for all n > Ni, Gj(X,n) > 0. If all the leading 
coefficients of all the Gj(X, ra)'s are "positive", there exists N' such that for all 
n> N', all the Gj(X, n)'s are positive. Therefore, Pi is nonterminating with input 
X' := A^'x. 

On the other hand, if Pi is nonterminating, does there exist an input X such 
that the leading coefficients of all the Gj (X, n) 's are "positive" ? We do not know 
the answer yet. However, with an assumption described below, the answer is yes. 

Assumption for the main algorithm: For any let" and any Cj k i(X, n), 
Cj k i2{X,n) being not identically zero implies 

mm(C jk w{X) + Cjkii(X) + C jk i 2 (X, y lu y 12 , ■ ■ ■ , y Sk i, y Sfc2 )) / 

subject to yl x + yf 2 = 1 (i = 1, . . . , s k ). 

It is not difficult to see that the assumption is equivalent to the following 
formula: 

[VXVy (C 3kl2 (X,n) = 0V 

Ai<i< Sfe 2/?i + vh = 1 - c jk io(x) + c jkn (x) + c jk i 2 (x, y) > o)] v 

{VX3Y (C jk i 2 {X,n) = 0V 

Ai<i< s , 2/?i + vl = 1 - C jk io(X) + C jkn (X) + C jkl2 (X, Y) < 0)] . 



Because Cj k i2(X,n) can be written as J2 fa(X) sin(nai27r) + fa(X) cos(nai27r), 



is/ 



where / is an index set, C jk i2{X, n) = is equivalent to /\ ig/ fn{X) = foiX) = 0. 
Thus, the assumption can be checked with real quantifier elimination techniques. 

Example 3. For those CV^'s in Example 2, let's check whether they satisfy the 
assumption. For clarity, we present a clear proof here rather than make use of 
any tool for real quantifier elimination. Take Gi(X, n) for example. It's clear that 
Ci2m{X, n) = if and only if xi = x 2 = 0. It is not difficult to compute 



D = inf Ci2o 2 (X,n) 



n>l 



x\ — X% + 6X1X2 \ ( lx\ — x\ — 2X1X2 



then 



If xi = X2 = does not hold, 5x i +x 2^ 2x i x 2 + jj < q because 

.5x1 +X% - 2X1X2,2 n 2 2 . 2 „ n2 . n 

( ) -D = — (5xi + x 2 - 2x1X2) < 0. 

Let 

M = min (Ci2oo(X) + Ci 20 i(X) + Ci 20 2(X, yu,yi 2 )), 

M = -(a;3+a;4) H h -D < 

if £1 = a; 2 = does not hold. Consequently, Gi(X, n) satisfies the assumption. 
Similarly it can be proved G 2 (X,n) and Gs(X,n) satisfy the assumption too. 
Thus the loop in Example 1 satisfies the assumption. 

We shall show in next section how hard it is to deal with the case that the 
assumption does not hold. Now, we are ready to describe our main algorithm. For 
the sake of brevity, the algorithm is described as a nondctcrministic algorithm. 
The basic idea is to guess a leading term for each Gj(X,n) first. Then, setting its 
coefficient be "positive" and the coefficients of the terms with higher order be 
"nonnegative" , we can get a semi-algebraic system (SAS). If one of our guess is 
satisfiable, i.e., one of the SASs has solutions, Pi is nonterminating. Otherwise, 
it is terminating. 

Algorithm Termination 

Step Compute the general expression of A n+m X. 

Step 1 Substitute A n+m X for X in P{X), and compute all Gj(X, n) (finite many, 
say, 3 = 1,...,L). 

Step 2 Guess a leading term for each Gj(X,n), say Cjh j i j r^ i r J k.. 
Step 3 Construct a semi-algebraic system S as follows. 

Sj = Cjkjij y A A(k>k j )\/(k=kj/\i>ij) Cjki(X,n) ^0, S = f\f =1 Sj. 



Step 4 If one of these systems is satisfiable, return "nonterminating" . Otherwise 
return "terminating". 



Remark 3. If the assumption for the main algorithm does not hold, then 
Termination is incomplete. That is, if it returns "nontcrminating" , the loop is 
nontcrminating. Otherwise, it tells nothing. 

Example 4- For the loop in Example 1, we have computed the Gt(X, n)'s in 
Example 2 and verified that it satisfies the assumption of the main algorithm 
in Example 3. We shall finish its termination decision in this example, following 
the steps in Termination. 

By Steps 2 and 3 of Termination, we should guess leading terms and con- 
struct SASs accordingly. To be concrete, let's take as an example one certain 
guess and suppose we have the following semi-algebraic system 

C C 120 (X,n) y 0, 
I C 2 2o{X,n) >- 0, 
[c 320 (X,n) >- 0. 

According to Notation [2] and Remark (5] the above inequalities are equivalent 

to (Vyn, 3/12) yii + yh = 1 — > 
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I ™ I 5i\<= I x 1 ^x 2 -vx 1 X2 . i x ± -x 2 -^1 x 2 n 
(£4 -t- -f) -\ | J J/ll i J J/12 > U. 

In Example 3, we have shown that 5x i+ x 2- 2x i x 2 -4- < o. Thus, the above pred- 
icate formula does not hold. In fact, none of the formulas obtained in Step 3 
holds. Thus, the loop in Example 1 is terminating. 

Remark 4- It is well known that real quantifier elimination is decidable from 
Tarski's work [15] . Therefore, the semi-algebraic systems in Step 3 can be solved. 
For the tools for solving semi-algebraic systems, please be referred to [617110118] . 

Remark 5. There are some techniques to decrease the amount of computation of 
the algorithm Termination. For example, we can use Lemma [3] when guessing 
leading terms for each Gj(X,n). 

Lemma 3. [4] Let £i, £2, • • • , £m £ C be a collection of distinct complex numbers 
such that = 1 and £i 7^ 1 for all i. Let oti, 02, • ■ ■ , a. m be any complex numbers 
and z„ = ai£i + . . . + a m Gk- Then one of the following is true: 

1. the real part Re(z n ) — for all n; or 

2. there is c < such that Re(z n ) < c for infinitely many n's. 

According to Lemma [3] if Cjkio = and Cj k i(X,n) is not identically zero 
w.r.t. n, then Cjki(X, n) can not be always nonnegative. According to the former 
discussion 

{G^ (X.n),...^-.! (X,n)} 

E Tfc+l £ T h +T 3 

k=l k=l 

are obtained from Pj (X, n) . Then the leading term with the greatest order among 
all the leading terms of the above set should not be of the form Cjkj i s n lj r£. where 



Cjk i o = 0. Thus, this should be avoided when guessing leading terms for each 

In the following subsections, we shall explain the details of the main steps of 
Termination and prove its correctness. 

3.3 Compute the Minimal Polynomials of a. + (3, a. — (3, a. • (3 and ^ 

At Step 1 of the main algorithm, we may compute rjk(l < k < M) which are the 
products of some £j's after substituting the general expression of A n+m X for X 
in P(X). In order to describe rjk, we need the minimal polynomials of r\h- In this 
subsection a method is presented to solve a more general problem. 

In [14] Strzebonski gave an algorithm to compute {a + (3, a — j3, a ■ (3, 
numerically where a and (3 are given algebraic numbers. Here we want to present 
another more intuitive method based on symbolic computation. 

In the following let a o (3 denote one of {a + /3, a — (3, a ■ /3, §}. Without 
loss of generality let's assume that the minimal polynomial of a is fi(x) whose 
degree is d\ and the minimal polynomial of (3 is f2(x) whose degree is di- We 
can bound a and (3 in W\ and W2, respectively, where W\, W2 are "boxes", by 
isolating the complex zeros of fi(x) and f2(x). Since the degree of a is d\ and 
the degree of (3 is c?2, the degree of ao/3 is at most dx-d^- Then there must exist 
X\, . . . x^-d.2+1 m 2 such that 

x x + x 2 a o (3 + . . . + x dl . d2+1 (a o /3) dl - d2 = 0. 

Thus we can design an algorithm to enumerate all the {x\, . . . ,Xd 1 -d 2 +i) €E 
jdi-di+i an( j chgck whether it is a solution. Since there exists one solution this 
algorithm must terminate and output a solution. Assume that its output is 
ao, 01, . . • , dd 1 -d 2 an d f(x) = aa + a\x + . . . + a dl .d 2 x dl ' d2 . Because f(a o (3) =0, 
the minimal polynomial of a o (3 is an irreducible factor of f(x). Factor f(x) in 
Q. Without loss of generality, we assume f(x) = gi(x) mi ■ g2(x) m2 . . . gd{x) md ■ 
We can check whether gj(x)(l < j < d) is the minimal polynomial of a o (3 by 
solving the following semi-algebraic system (SAS): 

{g 3 (xoy) = 0,/i(a:) = Q,f 2 (y) =0,xeWt,yG W 2 }. 

If it is satisfiable, gj{x) is the minimal polynomial of a o f3; otherwise it is not. 
Thus the minimal polynomial of a o [3 can be obtained. 

3.4 Check Whether the Argument of a Is a Rational Multiple of it 

At Step 1 of the main algorithm, T 3 is necessary for defining Gj(X,n). Thus for 
a given r/k we have to check whether its argument is a rational multiple of ir and 
if it is, we need to know the period of t^V This subsection aims at this problem. 

Suppose the minimal polynomial of a is p(x) whose degree is d. Without 
loss of generality suppose a — re 02nl . We can bound a in W by isolating all the 
complex roots of p{x). Since the degree of a is d, the degree of a must be d. 



Then the degree of a ■ a = r 2 is at most d 2 . Thus the degree of r is at most 2d 2 . 
The degree of a -1 is at most 2d 2 because the degree of a -1 is the same as the 
degree of a. Since the degree of a is d the degree of a ■ r^ 1 — e 4/327r is at most 
2d 3 . 

If (3 is a rational number, a must be a unit root and its minimal polynomial 
must be a cyclotomic polynomial. As a result if /3 is a rational number the mini- 
mal polynomial of a must be a cyclotomic polynomial whose degree is less than 
or equal to 2n 3 . All the cyclotomic polynomials can be computed explicitly ac- 
cording to the theory of cyclotomic field. Thus let CPj{x) denote the cyclotomic 
polynomial whose degree is j. Then, (3 is a rational number if and only if the 
following is satisfiable 

2d 3 

3r((r ji 0)/\(\/ CP 3 {x/r) = 0) f\(p(x) = 0) /\(x e W)). 

3=1 

Because CPj(x/r) = •<=>• r^CPj(x/r) — 0, the above quantifier formula is 
decidable. If it's satisfiable the minimal polynomial of a can be computed by 
checking whether 

3r((r ^ 0) (x/r) = 0) /\(p(x) = 0) e W)) 

is satisfiable one by one. 

3.5 Check Rational Independence 

Given a set of algebraic numbers, ai = e /5l2,rl , . . . ,a d = e^ d2vl , where /3i, . . . ,/3 d 
are irrational numbers. In this subsection we present a method to check whether 
... ,f3d are rationally independent. 

Definition 4. Irrational numbers . . . ,(3 d are rationally independent if there 

d 

does not exist rational numbers ai, . . . ,a d such that a jPj £ Q- 

3=1 

Obviously, /3i,...,/3d are rationally independent if and only if l,/3i, . . . , /3<j 
are linearly independent in Q. It can be deduced that /3i , . . . , (3d are rationally 

d 

independent if and only if V(6i, . . . ,b d ) e Z d , J2 bjPj $ Z. 

3=1 

Lemma 4. [T] Let Xi , . . . , A m w'f/i m > 2 be linearly dependent logarithms of 
algebraic numbers. Define ay = e Aj (l < j < m). For 1 < j < m, let log A, > 1 be 
an upper bound for max{ft(ay), ^y-} where D is the degree of the number field 
K = Q(ai, . . . ,« m ) over Q and h{a) denotes the absolute logarithmic height of a. 
Then there exist rational integers ni,... ,n m , not all of which are zero, such that 
m\i + ... + n m \ m = and\n k \ < (ll(m - l)!) 3 )" 1 ' 1 (^ A ^-^A m ) j orl < fc < m . 

Remark 6. Baker is the first one to use his transcendence arguments to establish 
such an estimate. However, the description here follows Lemma 7.19 in [17j . 



Sequence {f3i, . . . , f3 d , 1} are linearly independent in Q iff {/3i27ri, . . . , p d 2ni, 2ni} 
are linearly independent in Q. If ni/3i27ri + . . . + n d (3 d 2-K\ + n d+ i27ri = 0, then 
e" lft2d • • • e"" 13 " 2 ^ = 1. That is a™ 1 ■ ■ • a£ d = 1. According to Lemma H we can 
decide whether {/3i, . . . , f3 d } are rationally independent by enumerating n k from 

[ _ {lldD y(losA 1 )---(logA d+1 ) 3 d (logA 1 )---(logA <i+1 ) 

log A k log A k 

for fc = 1, . . . , d and checking whether a" 1 ■ ■ ■ a d d = 1. For any given (ni , . . . , n d ), 
whether a™ 1 • • ■ a 1 d >d = 1 can be determined by checking whether the following 
SAS has solutions. 

{x?—x2* = i, ©(xj) = o, xjeWj, j = i,...,d}, 



where qj is the minimal polynomial of ay and Wj contains only one complex 
root of qj for j = 1, . . . , d. 

3.6 Compute the Infimum of Cjki2 

To prove the correctness of our main algorithm, we need some further results. 

First, let's introduce a lemma in ergodic theory. Let 5* be the unit circumfer- 
ence. Usually any point on 5, say (a, 6), is denoted as a complex number a + bi. 
Define tt : R m -» S m as 

(xi,...,x m ) -> (e ,...,e m ). 
Define m-torus T m = {(e 112 "', . . . , e^ 2 ^ Xj <= R} and L w(o0 : T m -> T m as 

^ e yi2vri e i/„ l 27ri-j +"1 )2?ri + <* m )2?ri -j 



Lemma 5. [TT] //a £ R m , the translation L^^(X) is ergodic iff for all K g Z m , 
{K,a) ^ Z where (K,a) stands for the inner product of K and a. 

According to Lemma [5] we know that if ct\ is an irrational number, the clo- 
sure of {e" Ql27rl }„>i is the unit circumference. Also, if qi, . . . ,a m are rationally 
independent, L W ^(X) is ergodic. Thus the closure of {£"( q) (0)}„>i is T m . 

Lemma 6. // a k i, . . . , a k3k are rationally independent and Cjkn is of the form 

Cj k i 2 {X,sm(na kl 2n),cos(na kl 2TY), . . . ,sm(na k3k 2n),cos(na kSk 2n)) 

for a fixed X , then 

inf {C jk i2(X, sm(na kl 2n),cos(na kl 2n), . .., sin(na k3k 2n), cos(na kak 2tt))} 

n>l 

is equal to mm{Cjki2(X,xi,yi, . . . ,x Sk ,y Sk )} subject to {x 2 + y 2 = 1,1 < i < s k }. 

Proof. According to Lemma [5] for any (x\ , yi , . . . , x Sk , y Sk ) there exists a subse- 
quence, say {ni}i>i, such that 

lim (sin(niQfei27r),cos(nia^i27r), . . . , sm(ma kSk 2TY), cos(ma k s k 27v)) 

i— >+oo 

= (xi,yi,...,x Sk ,y Sk ). 



Theorem 3. Let 7, = (nTji +j )a ki 2ir (1 < i < s k ) and suppose that 
Cjki2 = Cjfci2(X,sin(7i),cos(7i),. . . , sin(7 Sfc ), cos(7 Sfc )). 

Then 

inf {Cjfei2(X,sin(7i),cos(7i),. . . , sin(7 Sfc ), cos(7 Sfc ))} 

n >1 

is equal to mm{Cj k i2(X, xi,yi,..., x Sk , y Sk )} subject to {x 2 + y 2 = 1, 1 < i < sk}. 

Proof. According to Lemma [5] and Lemma [HI it's sufficient to prove that T Sk is 
the closure of {(e 711 , . . . ,e 7s fc')} n >i. Because {a>i, . . . ,a 3k } are rationally indepen- 
dent, {Tjiai, . . . ,Tj/Q Sfc } are rationally independent, too. Thus, T Sk is the closure 
of 

{(e n7 >' a « 2,ri ,... J e n ' r >' a! W 2 ' ri )}„> 1 . The result of rotating 

(c nT 3' akl2lri e nT i' aks k 27Ti ) 

by (j «fti27r, j a kak 2n) is (e 71 ', . . . je 7 ^ 1 ). Consequently, T" k is the closure of 
{(e 71 ', . . . ,e 7s t ')} n >i. That completes the proof. 

3.7 Correctness 

For each j, Pj(X,n) can be written as 

D jl0 (X,n)r? + D jn (X,n)nr? + ... + D Jldl (X,n)n dl r? + ••• 
D 3 m)(X,n)r r H + Dj H i{X,n)nrH + . . . + D jH d H {X,n)n dH r^. 

The DjuS in the above are real because Pj(X,n) G R and those n l r k 's are of 
different orders. Just like Cjki, P>jki can be divided into three parts, 

D jk i = D mo {X) +D jkn (X,n) + D ]kl2 {X,n). 

Because D jM0 (X) contains no e (nT > +j " '^"'s, the e (- nT j'+i") a i a " l > s contained in 
Djkn(X, n) are periodic and the e (nT j' +J ) a i 2wlJ s contained in D jk i2(X, n) are not, 
Djkio, Djkn and Dj k n are all real. Since Cj k i(X,n) results from Djki(X,n), we get 
the following lemma. 

Lemma 7. For alln G N and each Cjki(X, n) = Cjkio(X) + Cjkii(X) + Cjkii(X, n), 
Cjkio(X) G E, C m i{X,n) G K and C* jfei 2(X,n) G K. 

Lemma 8. Consider C jk i{X,n) = Cjkio(X) + Cjkii(X) + C jk i2{X,n). 

1. inf {Cjki(X, n))} > if and only if 

mm{Cjkio(X) + Cjkii{X) + Cjkia{X,yu,yi2, . .., y 3k i, Vs k 2)} > 

subject to {y 2 ! + y 2 2 = 1, 1 < t < s k }. 

2. If I = Cjkio(X) + Cjkii(X) + D < 0, then there is c < such that C jk i{X,n) < c 
/or infinitely many n 's, where 

D = min{C jfe i2(X, yii,2/i2, . . . , y Bk x, ys h %)} 

subject to {j/ 2 i + y 2 t2 = 1, 1 < t < s k }. 



Proof. 1. It has been proved that 

inf Cj M2 (X,n) = min{C jk i 2 {X, y n , yn, . . ■ ,y Sk i,y Sk 2)} 

subject to {yti + y\ 2 = 1, 1 < t < s k }. Consequently 

inf {Cjki(X,n))} = min{C'. )fc ioPO + Cj k n(X) + G jk i 2 {X, yn, y 12 , . . . ,y Sk i,y Sk2 )}. 

rt> 1 

subject to {j/ t 2 ! + 2/f 2 = 1, 1 < t < s h }. 
2. Let Y = (2/11,2/12, •• • ,ys k i,y Sk 2) and V' = (yn,yi 2 , ■ ■ • ,2/^1,2/^2)- £> can be at- 
tained because {y n + yl 2 = 1, 1 < t < s k } is abounded closed set and C jk i 2 (X, Y) 
is a continuous function of Y. Assume that D = Cj k i 2 (X, Y'). Since I < there 
exists a neighborhood of Y', say U, such that 

VY £(/=>■ C 3fci0 (X) + C jfci i(X) + C jM2 (X, Y) < J/2. 
Let c = J/2 and 7, = (nTj, + j )a ki 2n. Because of the density of 

there are infinitely many n's such that (cos 71, sin 71, . . . , cos 7 S)t , sin 7 Sfe ) lies 
in U. Thus there are infinitely many n's such that Cj k t(X,n) < c. 

If the main algorithm, Termination, finds one solution, Xo, the leading coef- 
ficient of Gj(Xo), say Cj k i(Xo, n), satisfies Cj k i{X,ri) y 0. According to the defini- 
tion of y there exist cj > (j = 1,...,L) such that Cj k i{x,n) > cj for all n. Thus 
Pi is nonterminating. This means that if the algorithm outputs "nonterminat- 
ing", then Pi is nonterminating indeed. 

On the other hand, if the main algorithm outputs "terminating" , then for 
any {C jkj i j (X,n),j = 1, . . . , L} there is a subset V C {1, . . . , L} such that 

/\ C jkjlj (X,n)yO 

is not satisfiable subject to 

AjgvCj'Mj >" 0A /\(k>k j -)v(k=k j Ai>t j ) c Jki{X,n) y 

According to the assumption for the main algorithm, we get that with 
the above constraints 

V7 G V, m£C Jkj iAX,n)<0. 

n>l 

Thus by Lemma [51 Vj € V,Cj k -i(X,n) is identically zero or there are infinitely 
many n's and some c < such that Cj k i (X, n) < c. That means Pi is terminating. 
Therefore, we get the following theorem. 

Theorem 4. Under the assumption of the main algorithm, Termination re- 
turns "terminating" if and only ifP± is terminating. 



4 Conjecture 



In this section, we shall discuss the general case of Pi wherein our assumption 
for the main algorithm may not hold. 

Suppose p(X) = p(xn, X12, ■ ■ ■ , Xmi, Xm2) € Q[X], and one of the loop condi- 
tions is p{X) 2 > 0. From the discussion in Section 3, we know that if we substi- 
tute A n+m X for X in the conditions, there must be a polynomial q such that the 
condition becomes q(X, sin(nai27r), cos(nai27r), . . . , sin(na m 27r), cos(na m 2-7r)) 2 > 0. 
Because p is arbitrary, q can be arbitrary. Further, on, ...,a m can be made ratio- 
nally independent because A can be arbitrary. It's not hard to see that we can 
construct a program Q such that it is terminating if and only if 

Sg,a = {n : (j , (sin(no;i27r), cos(nai27r), . . . , sin(na m 27r), cos(na m 2-7r)) = 0} 

contains infinitely many elements. 

For any p(X) e Q[X] the decision problem "whether {X ■ p(X) = 0} f| % 2m = 0" 
is undecidable. Z 2m is a "regular" set while 

£ = {(sin(nai27r), cos(nai27r), . . . , sin(na m 27r), cos(na m 2-7r))} n >i 

is a chaotic set when {ai ,a m } are rationally independent according to the er- 
godic theory. Intuitively, deciding whether § Pt0l = {X : p(X) = 0} f]£ = is more 
difficult than deciding whether {X : p(X) = 0}f)Z 2m = 0. So, we intuitively 
guess the decision problem "whether S p , a is empty" is undecidable. Following 
the same idea, we guess the decision problem "whether S Pja contains infinitely 
many elements" is much more difficult and thus undecidable. Thus, we make the 
following conjecture: 

Conjecture. The decision problem "whether the loop Pi is terminating over 
R" is undecidable. 

5 Conclusion 

In this paper we have proved that termination of Pi over Z is undecidable. 
Then we give a relatively complete algorithm, with an assumption, to determine 
whether Pi is terminating over R. If the assumption holds, Pi is terminating 
iff our algorithm outputs "terminating" . If the assumption does not hold, Pi is 
nonterminating if the algorithm outputs "nonterminating" . We demonstrate the 
main steps of our algorithm by an example. Finally we show how hard it is to 
determine the termination of Pi by reducing its termination to the problem of 
" whether § p , a has infinite many elements". We conjecture the latter problem 
is undecidable. Thus, if our conjecture holds, the termination of Pi over R is 
undecidable. 

Acknowledgement 

The authors would like to thank Prof. Lu Yang for proposing the problem to 
us and encouraging us to work on it. We thank Prof. Michel Waldschmidt for 



his introduction to Baker's work. We are also grateful to Prof. Bob Caviness, 
Prof. Hoon Hong and Prof. Daniel Richardson for their insightful suggestions. 
Finally we would like to thank Prof. Chaochen Zhou, Prof. Naijun Zhan and all 
the other members in our group for the discussion at our weekly seminar. 

References 

1. A. Baker: Linear Forms in the Logarithms of Algebraic Numbers I, II, III, IV. 
Mathematika 13: 204-216, 1966; ibid., 14: 102-107, 220-228, 1967; ibid., 15: 204- 
216, 1968. 

2. V.D. Blondel, O. Bournez, P. Koiran, C.H. Papadimitriou and J.N. Tsitsiklis: De- 
ciding stability and mortality of piecewise affine dynamical systems. Theoretical 
Computer Science, 255(1-2): 687-696, 2001. 

3. A.R. Bradley, Z. Manna and H.B. Sipma: Termination of Polynomial Programs. In 
Proc. Verification, Model- Checking, and Abstract-Interpretation (VMCAI) January 
2005, LNCS 3385, pp. 113-129, 2005. 

4. M. Braverman: Termination of Integer Linear Programs. CAV 2006, LNCS 4114, 
pp.372-385, 2006. 

5. Y. Chen, B. Xia, L. Yang, N. Zhan and C. Zhou: Discovering Non-linear ranking 
functions by Solving Semi-algebraic Systems. LNCS 4711, pp.34-49, 2007. 

6. G. E. Collins: Quantifier elimination for real closed fields by cylindrical algebraic 
decomposition. In: Automata Theory and Formal Languages (Brakhage, H., ed.), 
LNCS 33, 134-165. Springer, Berlin Heidelberg, 1975. 

7. G. E. Collins and H. Hong: Partial cylindrical algebraic decomposition for quanti- 
fier elimination, Journal of Symbolic Computation, 12: 299-328, 1991. 

8. M. colon and H.B. Sipma: Synthesis of linear ranking functions. In TACAS01, 
LNCS 2031, pp.67-81, 2001. 

9. D. Dams, R. Gerth, and O. Grumberg: A heuristic for the automatic generation 
of ranking functions. In Workshop on Advances in Verification (WAVeOO), pp. 1-8, 
2000. 

10. A. Dolzman and T. Sturm: REDLOG: Computer algebra meets computer logic. 
ACM SIGSAM Bulletin, 31(2): 2-9, 1997. 

11. R. Mane: Ergodic Theory and Differentiable Dynamics. Springer- Verlag, NewYork, 
1987. 

12. A. Podelski and A. Rybalchenko: A complete method for the synthesis of linear 
ranking functions. In VMCAI, LNCS 2937, pp.465-486, 2004. 

13. R. Stanley: Enumerative Combinatorics, vol. 1. Cambridge University Press, 1997. 

14. A.W. Strzebonski: Computing in the Field of Complex Algebraic Numbers. J. 
Symbolic Computation 24: 647-656, 1997. 

15. A. Tarski: A Decision Method for Elementary Algebra and Geometry (2nd edn.). 
University of California Press, Berkeley, 1951. 

16. A. Tiwari: Termination of Linear Programs. In R. Alur and D.A. Peled(Eds.): CAV 
2004, LNCS 3114, pp.70-82, 2004. 

17. M. Waldschimidt: Diophantine Approximation on Linear Algebraic Groups. 
Springer- Verlag, Berlin, 2000. 

18. B. Xia: DISCOVERER: A tool for solving semi-algebraic systems, Software Demo 
at ISSAC2007, Waterloo, July 30, 2007. Also: ACM SIGSAM Bulletin, 41(3), 102- 
103, 2007. 



